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ABSTRACT 



Aims. Studies of the debris disc phenomenon have shown that most systems are analogous to the Edgeworth-Kuiper Belt (EKB). 
In this study we aim to determine how many of the IRAS 25/um excesses towards A stars, which may be indicative of asteroid belt 
analogues, are real, and investigate where the dust must lie and so build up a picture of what these systems are like. 
Methods. We observe using ground-based mid-infrared imaging with TIMMI2, VISIR, Michelle and TReCS a sample of A and B- 
type main sequence stars previously reported as having mid-infrared excess. We combine modelling of the emission spectrum from 
multi-wavelength photometry with a modelling technique designed to constrain the radial extent of emission in mid-infrared imaging 
to constrain the possible location of the debris. 

Results. We independently confirm the presence of warm dust around three of the candidates: HD 3003, HD 80950 and n Tel. For 
the binary HD3003 a stability analysis indicates the dust is either circumstellar and lying at ~ 4 AU with the binary orbiting at > 
14AU, or the dust lies in an unstable location; there is tentative evidence for temporal evolution of its excess emission on a ~20 
year timescale. For 7 of the targets we present quantitative limits on the location of dust around the star based on the unresolved 
imaging. We demonstrate that the disc around HD71 155 must have multiple spatially distinct components at 2 and 60AU. We model 
the limits of current instrumentation to resolve debris disc emission and show that most of the known A star debris discs which could 
be readily resolved at 18/jm on 8m instruments have been resolved, but identify several that could be resolved with deep (> 8 hours 
total) integrations (such as HD19356, HD139006 and HD102647). 

Conclusions. Limits from unresolved imaging can help distinguish between competing models of the disc emission, but resolved 
imaging is key to an unambiguous determination of the disc location. Modelling of the detection limits for extended emission can be 
useful for targeting future observational campaigns towards sources most likely to be resolved. MIRI on the JWST will be able to 
resolve the majority of the known A star debris disc population. METIS on the E-ELT will provide the opportunity to explore the hot 
disc population more thoroughly by detecting extended emission down to where calibration accuracy limits disc detection through 
photometry alone, reaching levels below 1 zodi for stars within lOpc. 

Key words, circumstellar matter - Techniques: high angular resolution - Infrared: stars 



1. Introduction 

Analysis of the IRAS database over the last 20 years has shown 
that over 300 main sequence stars have dust discs around them. 
This material is thought to be the debris left over at the end o f 
the planet formation process (e.g. iMannin gs & Barlow 1998). 
The spectral energy distribution (SED) of this excess in the best 
studied cases (e.g., Vega, f3 Pictoris, Fomalhaut, e Eridani) peaks 
longward of 60/mi implying that this dust is cool (<80K), and 
so resides in Edgeworth-Kuiper belt (EKB)-like regions in the 
systems. The EKB-like location and analogy is confirmed in the 
majority of cases w here these discs have been resolved (in th er- 
mal emission, e.g. iHolland et alj 19981: iGreaves et aL 20051; in 
scattered l ight imaging, e.g. iKalas et al .1 120071: Schneider et al. 
2006, l2Q09h . since the dust is shown to lie >40AU from the 
stars, and its short lifetime means that it must be continually 
replen ished by the cqllis ional destruction of km-sized planetes- 
imals ( Wvatt & Den1l2002h . The inner 40AU radius hole is thus 
thought to arise from clearing by an unseen planetary system, 
the existence of which is supported by the presence of clumps 
and asymmetries seen in the structure of the dust rings ( e.g. 
IWvatt et al.|[l999HGreaves et al.lll998l:rAugereau et al.ll2001l) . 



On closer examination the inner regions of known cool de- 
bris disc systems are much more complex than simply 'dust- 
free inner holes'. f3 Pictoris has (a r elatively small amount 
of) resolved dust in this inner regio n (Lagage & P antinl \l99^c 
iTelesco et afl 120051: iBoccaletti et all 120 09). thought to be there 
because this is a young (12Myr, Zu ckerman et al.l 1200 ll) tran- 
sitional system in which these regio ns have yet to be fully 
cleared by planet formation processes. lAbsil et al] (|2006) have 
recently presented interferome tric data showing Ve ga (thought 
to be around 380-500Myr old. iPeterson et al.ll2006h is likely to 
possess extended dust emission within 8AU, with evidence for 
a sim ilar w arm dust propulation around Fomalhaut (Absil et al. 
2009), and d i Folco et al.l (|2007) have also recently presented 
evidence for hot dust around the lOGyr old t Ceti. Multiple- 
component discs mirror our own Solar System, with our debris 
disc concentrated in the asteroid belt and EKB. It is possible that 
more sources with known cold discs have dust in the inner sys- 
tem, but the difficulty of separating hot dust emission fro m the 
stellar photosphere often limits detection. lChen et aU d2009l) pre- 
sented a sample of 1 1 debris systems believed to have multiple 
dust belts, although confirmation of the multiple components via 
resolved imaging is required for several of these stars. 
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The presence of excess in the mid-infrared range is of partic- 
ular interest. The temperature of dust emitting at 24yum suggests 
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that it should lie close to the star, in regions of a few to a few 
tens of AU. These are regions in which we might expect planets 
to reside (see above), and so the origin of the emission and pos- 
sible links with any cold debris in the system must be explored. 
Surveys suggest that 24yum excess may be common, with around 
half of main sequence stars t hat exhibit excess mid-inf rared 
emission in the IRAS database (Mannings & Barlow 1998) hav- 
ing an excess at 25fim only dZuckermanl 12001). However, ex- 
cesses takenjrom thelRAS database cannot be used at face 
value. ISong et aD d2002l) . who searched the IRAS database for 
excess emission towards M-type stars, noted that when search- 
ing a large number of stars for excesses close to the detection 
threshold, a number of false positives must be expected due to 
noise. There have also been a few instances in which the IRAS 
excess has been shown to be attributed to background objects 
that fall within the relatively large IRAS beams (>30"; such ob- 
jects range from highly reddened ca rbon stars or Class II YSO's 
iLisse et al1l2002l to distant galaxies fSheret et afll2004l) . Another 
possible so urce of mid-infra red excess emission is reflection 
nebulosity dKalas et al J 12002). Indeed it is now routine for pa- 
pers discussing the excess sources found by IRAS to address the 
possibility that some of these are bogus debris discs ( Moo r et aD 
2006; Rhe e et al.ll200l . 

Recent statistical studies using Spitzer data have revealed a 
large sample of A s tars with exces s emission at 24 and/or 70/itn 
dRieke et al.l l2005t ISu et al] 12006). Spitzer observations bene- 
fit from smaller beamsizes and higher resolution compared to 
IRAS, improving the reliability of excess measurements. The 24 
and 70//m excesses around A stars have a wide variation in levels 
amongst systems of similar ages, but overall there is a decrease 
in the upper envelope of excess i nversely proportional to time 
dRieke et al .1120051: ISu et al . 1 120061) . These features can be inter- 
preted in terms of a steady-state evolution of belt-like planetesi- 
mal discs in a collisional cascade, where the fall-off with time is 
due to the collisional grinding away of material and the variation 
in excess levels between systems of similar ages can be repro- 
duced by_^ariationsjn initial disc mass and plan etesimal belt 
radii (Wyatt et al. 2007b). Although recent work bv lCurrie et al.l 
(2008) has shown that excess emission at 24/mi around A stars 
increases from 5-10 Myr and peaks around 10-15 Myr before 
declining with age, which is not predicted by the steady-state 
model, this can be explai ned by the delayed formatio n of Pluto- 
sized bodies in the disc (Ke nvon & Bromlevll2004al) . It is only 
when Pluto-sized bodies are formed that the orbits of planetes- 
imals are stirred to high relative velocities and the steady-state 
collisional cascade can begin. 

However this model does not account for the possibil- 
ity of multiple belts (see above). There have also been other 
models proposed to explain the variance in excess levels be- 
tween similarly aged stars and the variety in the SEP slopes 
of tho se stars with excess at both 24 and 70//m in iRieke et ail 
d2005l) . which can be int erpreted as evidence for stochastic evo- 
lution dRhee et al.l l2007). A collision similar to the Earth-Moon 
forming massive collision has been proposed as the explana- 
tion for the large excess and spectral features observed around 
HD172555 dLisse et al.ll2009l) . Similarly short-lived dust pro- 
duction origins have been proposed for dust found close to the 
central star around several Sun-like stars. For these systems a 
steady-state collisional cascade production from a spatially co- 
incident planetesimal belt cannot exp lain the levels of excess 
dWvatt et al.ll2007at lLohne et al.l l2008). For debris discs with a 
cold dust population in addition to hot dust emitting at 24/im 
the outer planetesimal belt could be feedin g the hot dust pop- 
ulation (as has been proposed for r\ Corvi. IWvatt et al.l l2007a: 



Smit h et al.ll2008l l2009bl) . However, the mechanism that might 
be transporting the dust from cold outer regions to the hot in- 
ner locations observed at 24yum is as yet unclear (possibilities 
include dynamica l scattering by a migrating planet, see e.g. 
Gom es et al ] |2005HBooth"et" al. 2009). 

To tackle these issues regarding the origin of the 24//m dust 
emission, and in particular explore what its presence may reveal 
about planet formation and as yet undetected planetary popula- 
tions, this paper looks at the constraints on the true dust distri- 
bution around a sample of A stars with 24/vm excess. This can 
be assessed indirectly from SED fitting to multi-wavelength in- 
frared photometry, but uncertainties arising from degeneracies 
in dust model fitting and the possibility of multiple temperatures 
of dust mean that determining radial location from SED fitting 
alone is challenging (see discussions of individual sources in 
section 4). Resolved imaging provides more direct constraints 
on the dust location. 

This paper is structured as follows; In $2] the sample selec- 
tion is described. In |3] we describe the various observational 
and analysis techniques employed for the observations, with the 
results and discussion of individual sources presented in §|4] An 
extension modelling technique is used to explore which of the 
A star discs in the literature may be fruitful subjects of future 
imaging in §0 The implications of these and the observational 
results are discussed in §|6] Conclusions are in |7] 



2. The Sample 

The sample consists of A and B stars with IRAS published de- 
tections of excess emission at 12 and/or 25 yumQ A first-cut was 
applied to the list of all published detections to produce a final 
sample of 1 1 candidates (Table[T]l. This first-cut consisted of the 
analysis outlined below to determine if the excess identified by 
IRAS was likely to be real. 

For e ach star J, H, and K band fluxes were obtained from 
2MASS ( Skrut skie et al. | l2006l) and V and B magnitudes from 
Tycho2 dH0g et al.ll2.Q00t) . The photosphe ric emission was then 
determined by adopting a Kurucz model dKuruczll 19791) for the 
appropriate spectral type (as listed in the Michigan Spectral 
Catalogues or SIMBAD) scaled to the K band flux. IRAS 
fluxes were extracted using SCANPI (the Scan Processing and 
Integration toolfl The expected stellar flux was multiplied by 
the colour-correction factor (at the levels described in the IRAS 
Explanatory Supplement) before subtraction from the IRAS 
flux to determine the excess. The proximity of the IRAS sources 
to the stars was also checked given the quoted uncertainty er- 
ror ellipse, since some surveys allowed excess sources to be 
up to 60 "offset and have si nce been shown to not be related 
(ISvlvester & Manningsll2000l) . 



3. Observations and Data Reduction 

3.1. Observations 

Observations were performed using: TIMMI2 on the ESO 3.6m 
telescope at La Silla (proposals 71. C-03 12, 72.C-0041 and 74.C- 
0070); VISIR on the ESO VLT (proposal 076.C-0305); and 



1 The sample stars are listed in the Debris Disc Database at 
http://www.roe.ac.uk/ukatc/research/topics/dust 

2 http://scanpi.ipac. caltech.edu:9000/ 

3 The IRAS Explanatory Supplement is available at 
http://irsa.ipac.caltech.edu/IRASdocs/exp.sup/ 
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Star name 


Stellar type 


Age 


Distance" 


IRAS fluxes (mJy)* 






Myr 


pc 


12/jm 


25pm 


HD 3003 


AOV 


50 c 


47 


148 + 16 (±20) 


34 + 241 (±20) 


HD 23281 


A6V 


626'' 


43 


251 + 36 (±20) 


59 + 107 (±28) 


HD 23432 


B8V 


100 e 


119 


ill + 230 (±43) 


29 + 1 1 5y (±42) 




AU V 


1 nnc 


in 


336 + 238 (±65) 


78 - 110 (±121) 


HD 38206 


AOV 


9' 


69 


130 + 7 (±23) 


30 + 67 (±21) 


A Gem s 


A3V 


560' 


29 


1166 + 214 (±28) 


271 + 161 (±52) 


HD 71155 


AOV 


169 c 


38 


675 + 162 (±31) 


157 + 249 (±41) 


HD 75416 


B8V 


5' 


97 


134 + 86 (±22) 


31 + 86 (±38) 


HD 80950 


AOV 


w 


81 


120 - 5 (±26) 


28 + 101 (±17) 


HD 141795 


Am 


450 E 


22 


1154 + 67 (±25) 


269 + 141 (±28) 


t] Tel' 


AOV 


12'' 


48 


263 + 138 (±34) 


61 + 394 (±22) 



Note s: "^Distance from pa rallax in Hipparcos; ''=Fluxes are shown as star + excess (±err or); c =Age taken from Song et al] J200lh : rf =Age taken 
fromlKunzli & North ( 1998); e = Age determined from theoretical evolutionary tracks of So ng et all d200 lh. see dis cussion in text; ^=Age taken 
from lRieke et al] d2005h : s HD 56537; ''= =Age from /? Pictoris association membership, Zuckerman et al. (2001); ' H D 181296. For those 
sou rces w ithout ages listed in the literature (HD 23432, HD 31295 and HD 141795) the p and C\ magnitudes as listed in[Hauck & Mermilliod 
d!998h were used to determ ine T^g and log(g) from the grids of M oon & Dwo retskv ( 1985). Then the theoretical evolutionary tracks of 
ISongetni ( 120011) were used to determine the approximate ages for these stars. 



Michelle and TReCS on the twin telescopes of Gemini (GN- 
2005B-Q-15 and GS-2005B-Q-67). All observations used chop- 
nod pattern to remove sky and telescope emission. A chop of 
10"in the North-South direction with a perpendicular 10"nod 
was used for the ESO observations. The Gemini observations 
used a 15"chop and parallel nod (also of 15") at 30° East of 
North. Observations of A Gem were performed with a chop at 
268° to ensure the image of the binary companion would fall on 
the array. 

This chop-nod pattern means that a simple co-addition of 
the data produces an image with two positive and two negative 
images of the source for the ESO observations, and one positive 
and two negative images at half the intensity level for the Gemini 
observations. A dark current offset is determined from median 
values for each row and column of the image (excluding pix- 
els on which the source image fell) and subtracted from the final 
frame. Pixels showing high levels of variation throughout the ob- 
servation (10 times the average) were masked off. Pixels show- 
ing very high or low gain (determined by comparing average sky 
emission detected across the image to that detected in each indi- 
vidual pixel detection) were also masked. In total an average of 
around 7% of pixels were removed in the TIMMI2 observations, 
and around 4% of pixels in the MICHELLE, TReCS and VISIR 
observations. Calibration observations of standard stars within a 
few degrees of the science object were taken immediately before 
and after science observations. The stan dards were chosen from 
the list of K and M giants identified bv lCohen etal.1 dl999h . In 
addition to photometric calibration, these standards were used to 
characterise the PSF and used for comparison with the science 
sources to detect any extension (see section [3~3"b . 

3.2. Photometry and Background/Companion Objects 

The multiple images resulting from the chop-nod pattern were 
co-added to get a final image by first determining the centroid of 
each of the individual images. Photometry was then performed 
using a 1" radius aperture for the TIMMI2 images and a CC'5 
radius aperture for the VISIR and MICHELLE images. These 
sizes were chosen to just exceed the full-width at half-maximum 
(FWHM) found for each instrument (average FWHM: 0'.'80 + 
0712 in the N band, and 1734 + 0710 at Q for TIMMI2; 07465 
+ 07161 at N on VISIR, and in the Q band a/597 ± 07166; for 



Michelle a/557 ± a/ 107 in N and a/579 + 07 101 in Q; and 
for the N band observations with TReCS a/475 + a/054). Note 
that the filters used in these observations were narrow band and 
so no colour-correction was applied. Residual statistical image 
noise was calculated using an annulus centred on the star with 
inner radius matching the outer radius of the aperture used for 
the photometry, and outer radius of twice the inner radius (so 2" 
for TIMMI2 and 1" for VISIR and MICHELLE). Typical levels 
for statistical noise at the 1 <x level in a half hour observation 
were 44mJy total in the 170 radius aperture of TIMMI2, 4 mJy 
and 12mJy for the a/5 aperture of VISIR in N and Q respec- 
tively, 6mJy in the a.' 5 aperture of MICHELLE and 3 mJy in the 
a/5 aperture of TReCS. Calibration uncertainty was determined 
from variation in standard star photometry, and was added in 
quadrature to statistical uncertainty to give the total error on the 
photometry as listed in Table [2] 



Smaller apertures were used to search for background 
sources and to place limits on undetected sources. The aperture 
radius was determined through examination of the standard star 
images. Circular apertures of increasing radius were centered on 
the standard star images and the radius giving the highest signal- 
to-noise (statistical noise only, as determined in the annuli listed 
above) was recorded. For each instrument and observing wave- 
length the median optimal radius for maximising signal-to-noise 
on the standard stars was chosen to optimise point source detec- 
tion. The apertures had radii of: a/8 for TIMMI2 observations; 
0/4 for MICHELLE and for TReCS; and a/32 and a/35 for the 
N and Q filters for VISIR. Apertures were systematically centred 
on each pixel of each array to search for > 3<x detections; where 
none were found the limits on any background object were based 
on the 3cr uncertainty in the aperture plus calibration uncertainty. 
For the non-photometric nights, limits were based on calibration 
to the IRAS flux of the object. The upper limits to background 
sources are listed in Table [2] These have been translated into a 
limit on the spectral type of any companion source. These spec- 
tral limits assume any companion is a main-sequence star at the 
same distance as the source, and are given as the hottest star that 
does not exceed the point source limits found in the imaging. 
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3.3. Extension testing and limits on disc size 

Evidence for extended emission was checked for all science tar- 
gets. The source's surface brightness profile was determined by 
calculating the average surface brightness in a series of annuli 
centred on the source of 2 pixel thickness by increasing inner 
radius from to 3", and this was compared to profiles of the 
standard targets. Finally the images of the point-like standard 
stars scaled to the peak of the science images were subtracted 
from the science images and the residuals checked for consis- 
tency with noise levels measured on the pre-subtraction image. 
A range of regions optimised for different disc geometries were 
tested for evidence of residual flux indicating spatially extended 
emission. These optimal regions were determined from exten- 
sive modell i ng wo rk and are outlined in detail in section 4 of 
ISmithetal.1 (120081) . 

To test the limits we can place on disc extension with un- 
resolved images, we performed the same PSF subtraction and 
residuals testing on models of stars + discs. The stellar compo- 
nent was modelled as a point source with flux set as predicted 
for each star in the appropriate filter (see Table Discs with 
different radii, width and inclination to the line of sight and with 
different levels of flux were added to the point source. The whole 
model was convolved with the PSF as modelled by standard star 
images to create a range of model images. Different standard star 
images were used to model the effects of PSF vari ation. This 
process is described in detail in Smi th et al.l (|2008). The disc 
geometries considered are simple ring-like discs with uniform 
brightness, with central radius r and width dr. The limits shown 
in Figures QJ4] are for disc widths dr = 0.2r (so a disc extending 
from r - 0.1 r to r + 0.1 r in radius) and dr - 2r (a disc from 
the central star to 2r). The central disc radius was varied from 
to 1'.'5 for the observations with VISIR and Gemini, and up 
to 4"for TIMMI2 observations. The flux of the disc was scaled 
from to 100% of the total flux of the source in the observed 
filter. Each model image was subjected to the same testing pro- 
cedures as the science image itself, testing regions of the point- 
source subtracted image that had been optimised for the detec- 
tion of extension for the disc parameters used as input. These 
optimised r e gions were based on modelling work described in 
ISmith et al.l 00081) . Re gions above the lines for different disc 
geometries in Figures Q]|4] represent disc models that were de- 
tected as extended objects (emission detected in optimal test- 
ing regions) at a level of at least 3<x (noise included pixel-pixel 
background noise and noise from PSF uncertainty as detailed in 
ISmith et al.l l2008). Regions below the lines (shaded area) repre- 
sent disc models that were not detected. The resulting limits on 
detecting extended emission are dependent on PSF stability for 
discs at small radii, and on the sensitivity of the observation for 
discs at larger radii (see Figures 2 and 3 of lSmith et al.ll2008T) . 

As we do not resolve any extended emission in the obser- 
vations presented in this paper, these limits are compared with 
the predicted radial location of the disc. The excess emission 
SED is fitted with a single-temperature blackbody which is con- 
verted to a radial offset assuming blackbody-like grains. If the 
grains which dominate the emission are small they will be in- 
efficient radiators, hotter than blackbody grains at a fixed radial 
location. Thus the predicted radial location will be an underesti- 
mate of the disc offset if small grains dominate the emission (as 
has been seen in scat tered light and thermal imaging of resolved 
discs, see section 5). Schneid er et al.l J2006) showed that for the 
HD181327 system its disc was imaged at a radius correspond- 
ing to 3 times that expected from a blackbody fit to the emission 
spectrum, a fact attributed to the emission from this disc being 



dominated by small inefficiently emitting but efficiently absorb- 
ing grains. Other resolved discs have been shown to have radii 
that can differ from the blackbody fit by up to a factor of 3 (see 
Table HJ. As the grain properties of the discs are unknown, in 
the results section the extension limits are compared to the ra- 
dius suggested by assuming blackbody grains and up to 3 x the 
blackbody radius. In all cases the extension limits are consistent 
with a disc lying at the blackbody radius (see Table [3} but the 
discs could also be dominated by smaller grains at a larger ra- 
dial offset (up to the extension limit). The exception in the single 
disc case is HD71155, for which the extension limits indicate 
the disc must have multiple belts (see discussion in section 4.2). 
We therefore also determine the minimum grain size that will 
not exceed the temperature fit to the excess emission when at the 
extension limit listed in Table|3] This calculation requires the as- 
sumption of grain composition which we take for reference to 
be non-porous grains with no ice inclusions with a silicate frac- 
tion o f 1/3 and 2/3 organic refractory material (see Wyatt & Dent! 
2002 for details of how grain temperatures were calculated for 
the assumed composition). These minimal grain sizes are listed 
as Z) gla i n in Table [3] In all cases these grains are smaller than the 
blowout limit (1.3^m for an AOV-type star assuming the grain 
properties given above, or greater for cooler stars), and so it's 
more likely that the true disc radius is smaller than the extension 
limit given in Table [3] 

4. Results 

We split our debris disc targets into 3 subgroups based on the ob- 
servational results: those with excess independently confirmed 
in our photometry; those for which limits can be placed on the 
extent of the disc with these observations; and those for which 
we cannot place limits on the disc with our data (HD 31295 and 
HD 38206). These categories are identified in Table [3] In ad- 
dition HD 23432 was found to have excess due to a reflection 
nebula and not a disc. This source together with HD 31295 and 
HD 38206 are described in Appendix A. The photometrically 
confirmed sources and those for which we can place limits on 
the disc with our observations are described below. 

4.1. Photometrically confirmed discs 

rj Tel: The excess emission towards rj Tel was confirmed in 
our N band photometry with TIMMI2 at a Act level of signif- 
icance (Table [2}. Resolved imaging of this target, revealing a 
two-component disc system with an outer component lying at 
24 AU resolved at 18yum and a further unresolved inner compo- 
nent which SEP fitt ing shows is at 4AU, is presented in detail in 
ISmith et aU d2009a). We shall not discuss this source further in 
this section, but will include this source in the discussion of the 
sample in section 6. 

HD 3003: This star was identified as having s ignificant 25/im 
excess (see Table [TJ by Oudmaiier et al. (1992). This star is a 
binary: both components are A stars with similar lu minosities 
(L A = 13.1L e ,L B = lO.OLM. lDommanget & Nvslll994l) . The B 
component was listed as being at an offset of 071, 143° East of 
North in 1925. The last confirmed observation of the separate 
componen ts was in 1964 wit h the B component at an offset of 
071, 171° dMason etal.ll200l . 

The star was observed with TIMMI2 at N and Q with follow- 
up on 8m telescopes in both bands (see Table|2]l. None of the ob- 
servations resolved the separate stellar components. Excess was 
confirmed photometrically at Q (detected 244+27mJy, expected 
113+2mJy from photosphere giving an excess of 131±27mJy 
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Table 2. The Observations 







Observation 




Predicted 






Results" 






Star 






Integration 




stellar 


Flux, 


Total 


Stats. 


Background limit 


name 


Filter 


/I, jj m 


time, s 


Instrument 


flux, mJy* 


mJy 


err., mJy 


err., mJy 


mly* 


Spc. lype' 


HD 3U03 


JNV.S 


9.56 


1800 


T| \ ir-\ jt I -> 

1 IMMlz 


426 


A O 1 

481 


5/ 


1 / 


S 40 


YSIYD.J V 




N12.9 


ion 

1Z.Z1 




TIMMI2 


262 


436 


64 


10 


<r 11 
S 2 / 






Q2 


1 S 7S 


^ 1 on 

J 1UU 


1 IMMlz 


1 12- 


oaa 
2Uu 


14U 


A A 

44 


< 148 






M-3 


11.66 


3400 




Zoo 


J /3 


34 


A 

4 


S 1U 








1 S 79 
i o. / z 


z^+uu 


V lolls. 






11 


q 

y 


S 21 




HD 2328 1 


N12.9 


19 91 
1Z.Z1 


19 no 


TT1V if A l\ 1 ~) 

1 IMMlz 


243 


339 


57 


10 


< 2.1 


Mjv 




SiC 


1 1 

1 l.OJ 




V1MK 


23 5 


232 


13 


2 


^ A 

< 4 






Q2 


1 C T9 
1 O. /Z 




VISIR 


104 


129 


9 


6 


< 13 




T TT~"v 11 ^11 

HD 23432 


N7.9 


"7 T7 


1 OAfl 

loUU 


T1MMI2 


299 


452 


192 


192 


^ ^ 1 
S j3 1 


AO V 




Nl 1.9 


1 1.J7 


1 son 

1 OUU 


1 IMMlz 


1 OA 

136 


141 


62 


14 


^ A A 

< 44 




HD 3 1 Zyj 


XT 1 

IN 1 


O.OU 


loUU 


1 IMMlz 


A /I Q 
045 


53 I 


in 
1 1 


13 


< 33 


M3V 


HD 38206 


N9.8 


9.56 


1800 


TIMMI2 


204 


196 


48 


16 


< 44 


G3.5V 


A Gem 


Nl 


8.60 


1800 


TIMMI2 


2244 


1997 


167 


22 


< 3D 


Mjjv 


Chop PA 


SiC 


1 1.5 J 


J4U 


VISIR 


1195 


701 


9 


4 


S / 




268° 


Q2 


1 £ 79 
10, /z 


1 sen 
1 00U 


VISIR 


484 


590 


71 





< 14 




Binary 


IN 1 


o.oU 


1 onn 
loUU 


1 IMMlz 


1 n 1 
1U1 


AO 

05 



5 








Chop PA 


SiC 


11. OJ 




VISIR 


54 


46 


4 








268° 


Q2 


1 S 79 
1 0. / z 


1 son 

1 OUU 


VISIR 


22 


31 


17 








HD /lljj 


JNz 


1U.05 


JUOU 


1 IMMlz 


53U 


12 so 


213 


Jo 


S 70 


JVH-. J V 




Si-5 


1 1 .ou 


1H-1U 


A /IT/^TTrT T 1 — ' 

MICHELLE 


722 


1050 


115 


12 


< 32 








1 c 

Io.jU 


9 1 (in 
zlUU 


1V11L- 11 EA-A-,c J 


OSA 
2o0 


Jyo 


QQ 

yy 




< 38 






Q2 


18.72 


3600 


VISIR 


280 


380 


36 


7 


< 16 






N1 1 

imi.7 


11.59 


1800 


TTMMT9 

1 IIVIIVIIZ 


1 AA 




00 


1 j 


<48 


F0V 


HD 80950 


N11.9 


11.59 


3600 


TIMMI2 


129 


175 


16 


10 


<25 


K4.5V 




SiC 


11.85 


1800 


VISIR 


124 


120 


23 


3 


< 7 






Q2 


18.72 


3760 


VISIR 


50 


119 


11 


9 


<20 




HD 141795 


N12.9 


12.21 


1200 


TIMMI2 


1115 


1138 


119 


19 


<49 


M4.5V 




Qa 


18.50 


1900 


MICHELLE 


491 


511 


39 


3 


< 8 




77 Tel 


N12.9 


12.21 


2500 


TIMMI2 


254 


351 


23 


10 


<25 


K2V 




Q2 


18.75 


1800 


TIMMI2 


109 


391 


182 


43 


< 136 





* The expected photospheric emission is determined by a Kurucz model profile appropriate to the spectral type of the star and scaled to the K 
band 2MASS magnitude as outlined in section 2 unless otherwise stated in the individual source description. Errors are lcr. M band TIMMI2 
observations were largely non-photometric and primarily used to improve pointing accuracy and thus are not listed in this table. Notes: " Errors 

are total errors (inclusive of calibration uncertainty and image noise). * Limits are 3cr upper limit to undetected object including calibration 
errors, or scaled to IRAS fluxes when conditions were non-photospheric. These limits are valid to within 28" of the detected source for TIMMI2 
observations, 12'.'6 for MICHELLE, 8'.' 1 for TReCS, and 1 Y'A of the source for VISIR observations. 1 The hottest spectral type of main-sequence 
star equidistant with the target that would have remained undetected within our field of view. See section 3.2 for details of this limit. 



- uncertainty on photosphere taken from fitting Kurucz model 
profiles to K band flux+lcr error, and uncertainty on the de- 
tected emission and on the photosphere were added in quadra- 
ture to give error on excess), although calibration uncertainty 
prevents confir mation at shorter w avelengths. This is in good 
agreement with ISmith et ail d2006l) who presented MIPS 24/vm 
confirmation of excess for this source (detected 155 + 9 mJy af- 
ter subtraction of the photosphere). The recent 24yum and 18/itn 
detections presented here suggest a lower excess than that sug- 
gested by IRAS photometry (275 + 20 mJy at 25 fim). Though 
the results are different at the Act level only, it is possible that 
the emission from the source may show temporal variance. A 
dust temperature fit of 265^K is consistent with the observed 
excess and shorter wavelength limits (see Table [3] and SED in 
Figure [TJ. Uncertainty on the blackbody temperature fit was 
taken from temperature range that allows a fit to within 3cr for 
all excess measurements. Undetected background targets within 
the TReCS field-of-view (section 3) are constrained to <10mJy 
at N, confirming the excess is centred on the star. 

The final images showed no evidence of extended emission 
at any band. The Q band VISIR images place the tightest con- 
straints on the disc size, putti ng a limit of <6.5AU on the disc's 
extent. If the binary listed in lMason~e t al. (200]]) is a true bi- 



nary at a separation of Cf.' l (4.7AU at 47pc, T abled) the stability 
analysis carried out by iHolman & Wiegertl (1 19991) would sug- 
gest the disc cannot be circumbinary, as the disc should be at a 
radius of > 0'.'24 (11.3AU or greater if the orbital separation is 
larger or the orbit is eccentric)0 We conclude that the disc must 
therefore be circumstellar. Assuming that the dust is around the 
primary star the temperature fit of 265K translates to an offset of 
4AU assuming black body grains using r = (278. 3 IT) 2 V^* (e.g. 
Back man & Parescd H993. where r is the dust location in AU, 
T is the temperture in Kelvin and L* is the stellar luminosity 
in units of L ). This assumes the grains are in thermal equilib- 
rium with their environment. Grains smaller than the wavelength 
at which the excess emission peaks are inefficient emitters and 
are thus hotter than blackbody temperature at that radial offset 



We have tacitly assumed here that the dust grains are distributed 
following the orbits of the parent planetesimals. As dust grains are af- 
fected by Poynting-Robertson drag they may have a different spatial 
distribution t o the parent bodies. However, following the equations in 
Wvatt (2005) we find that the grain collisional lifetime (f co n) is shorter 
than the Poynting-Robertson drag timescale (/ PR ) for all disc radii up 
to the resolution limit, i.e. that f co u < 0.04/pr for r < 6.5AU. The dust 
grains are therefore likely to be collisionally dominated and occupy a 
spatial distribution similar to the parent population. 
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Table 3. The Results 
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90 
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1.59 
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0.012 


5.64 


0.002 


HD 23281 


210 


5.36 


0.12 


< 9 5+^ 


0.68 


3.82 


0.024 


HD 75416 


250 


11.1 


0.11 


< 55.2 + - 




7.98 


2.44 


HD 141795 


250 


4.64 


0.22 


< 6 l +0 - 3 


0.73 


4.43 


0.008 


Unconstrained by these observations 










HD 31295 


80 


52.6 


1.4 






5.85 


9.216 


HD 38206 


90 


48.4 


0.70 






15.1 


84.96 



Note that the objects with no extension limits have too low a fractional excess for the extension to have been detected in the images regardless of 
size. Estimates of radius are based on blackbody fits and could be up to three times larger if the grains responsible for the emission are small. Limits 
on extension shown here are for a narrow face-on discs. Errors on this limit arise from 3 cr photometric errors and errors on the determination of 
the photospheric emission at the observed wavelength. These sources are combined to give 3<x uncertainties on F disc , which in turn give a range 
of radii at which a disc could have been detected as extended emission. Horizontal lines indicate division into photometrically confirmed debris 
discs, and sources for which our results provide constraints on the discs (sections 4. 1 and 4.2 respectively). 

Notes: " See section|6]for details of this limit; b the minimum grain diameter (for the composition assumed in section 3.3) such that grains at the 
maximum radius given by the extension testing limits do not exceed the temperature of the blackbody fit to the excess emission (where limits on 
the disc radius are large compared to the blackbody radius no limit can be placed on the mini mum grain size, as e ven the smallest grains cannot 
reach the fitted temperature at the extension limit radius); c These disc locations are taken from Smith et al. (2009a); d binary detected. 



from the star. Such grains can be offset by three times the ra- 
dius suggested by a blackbody approximation (see discussion in 
section 3.3). Assuming a radius of 4 AU there are then two pos- 
sibilities for the system: either the dust is in a stable location and 
the binary must have a semi-major axis of at least 14.4AU (or 
larger if the binary orbit is eccentric, according to equations of 
iHolman & Wiegertlll999h : or the binary is closer to the star and 
the dust is unstable. Such unstable dust population s have already 
been detected in a small number of binary systems (Trillin g et alJ 
120071) . However, if the dust is in an unstable region this could 
naturally explain the tentative evidence for temporal evolution 
in the level of excess. If the dust resides in a stable location, 
the motion of the binary changing the overall illumination of the 
system as it travels on its orbit could also possibly explain any 
temporal evolution. Determining the orbit of the binary will be a 
crucial step in determining the stability of the dust in this system. 

HD 80950: This star was identified by Mannings & B arlowl 
(1998) as a possible host of mid-infrared excess based on the 
IRAS 25yum measurement of its flux (excess 101 + 17 mJy, see 
TableQ}- The source was observed with TIMMI2 at N and VISIR 
at N and Q, with the Q band photometry allowing a confirma- 
tion of the excess (total detected flux 119+11 mJy, predicted 
stellar photospheric emission in this filter 50±2mJy with uncer- 
tainty taken from 2MASS Kband uncertainty, see Table 0. A 
fit to the Q band ex c ess em ission and the 24 and 70/mi excess 
reported in ISu et at] d2006l) suggests a temperature of 180*?qK 
for the excess emission (uncertainty determined by range of 
temperature fi t s that fit all excess measurements to within 3<x). 
Mora les et all (2009) used a similar blackbody temperature of 
188K to fit the Spitzer MIPS and IRS data on this target. No 
background/companion sources were detected in any of the im- 
ages (see Table|2]for brightness limits on such sources). No ev- 
idence of extension was detected in any of the images. The re- 



sulting limits on possible disc sizes and geometry are relatively 
broad due to high levels of variation in the PSF during the ob- 
servations. For face-on discs the limits suggest a disc radius of 
<24.5AU if the dust is distributed in a narrow ring (at a distance 
of 81pc, TableQ}, or <61AU for abroad face-on disc (see Figure 
[TJ. Assuming the grains are blackbody-like suggests an offset 
of 13.6AU, consistent with these limits for all disc geometries. 
Grains much smaller than the peak of emission can lying at 3x 
blackbody radius (see discussion of HD 3003 and section 3.3) 
would lie at 40.8AU (075 1) which should be detectable in 8m ob- 
servations with a more stable PSF. With our current constraints 
we can say that if very small grains dominate the emission, the 
dust must be distributed over a broad radial range. 

4.2. Limits on disc extension 

HD 71155: This star was identified bv lCotel (1 19871) as being a 
host of mid-infrared excess based on IRAS observations (see 
Table [TJ- Calibration uncertainty prevented photometric confir- 
mation of the excess in our observations, but we can rule out 
background objects within the TIMMI2 field of view of >98mJy 
at N or within the VISIR field-of-view >16mJy at O (Table 13 
We thu s co nfirm any excess should lie on the target. Riek e et al.l 
d2005l) and ISu et alJ (120061) used MIPS photometry to confirm 
the excess at 24 and 70/im. These results suggested somewhat 
lower excess than found with IRAS, but the difference is not 
significant at 3<x. In the following analysis we retain the MIPS 
24/vm result (detected 302±19mJy) in preference to the IRAS 
25yum result to take advantage of the reduced errors. 

The more recent measurements (upper limits at N and Q pre- 
sented in this paper and the MIPS excesses) allow a single tem- 
perature dust fit (model A) at 120+15K (uncertainty from all 
temperatures that fits excess data with 3<x, see Tableland Figure 
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Fig. 1. The observational results for sources with excess emission photometrically confirmed by our observations. Top: The results 
for HD 3003. The SED (left) includes a fit to the excess emission at a temperature of 265K (dotted line). The dashed lines indicate 
the photospheres of the two A star components, with the solid line being the sum of their emission. Bottom: The results for HD 
80950. The excess emission is fit with a temperature of 180K (dotted line - solid line shows photospheric emission). For both targets 
the limits on disc extension are shown in the right-hand plots, with the shaded region indicating the area the disc must lie in to have 
avoided being resolved. The difference between the limits on disc extension arise from the differences in the PSF model for the 
different observations. Specifically, the standard star observations associated with HD80950 showed large changes in the wings of 
the PSF over the course of the observations, giving rise to large uncertainties in the detection of excess emission extended over a 
broad spatial region. Thus the limits on broad disc structures (dashed lines) with small central radius are poorer for HD80950 than 
for HD3003. The disc flux marked on these plots comes directly from the VISIR photometry. The radial location marked arises 
from assuming the emitting material is blackbody-like. These results are discussed in detail in section 4.1. 



2 tem perature similar to the 105K fit suggested by ISu et al.l 
2006). However, if the 12 fim excess detected by IRAS is taken 
into account, which is not ruled out by the limits presented here, 
then a two-temperature dust model (model B) fits the spectral 
energy distribution better (500+180K and 90+20K, see Figure 
. Although no evidence for extension was seen on any of our 
images, the non-detection in the VISIR Q band image allows 
us to place constraints on the dust location in the context of 
the two alternative models making the assumption that the dif- 
ferent temperatures represent different radial locations. Figure|2] 
shows that we should have detected resolved emission if the dust 
was located at a radius of 34AU (assuming blackbody grains at 
120K; model A), regardless of disc geometry. The results there- 
fore sup_2or^Jwo_£omponent disc model, as was found for 77 
Tel (ISmith et al.l l2009a). The inner component of this model is 
limited to <8.4AU assuming a face-on orientation (Figure 0. 
The predicted location of 500K dust assuming blackbody-like 
gr ains is ~2AU. This pre dicted location agrees with the results 
of M oerchen et al.l d2009l) . who found extended emission around 
this source at 10.4/im consistent with a disc at 2AU. The outer 



component is predicted to be at 61AU (1'.'59) where sufficiently 
deep observations on current instruments could resolve this disc. 

A Gem: This source was listed in ICheng et al.l d 1992b . a 
study of main-sequence A-stars, as having an IRAS excess. This 
source is notable as having one of the largest 24/mi excesses 
amongst older stars dRieke et al.l 12005; stellar age is 560Myr). 
The star is listed in the Washington Double Star Catalogue as 
having a visual binary companion at a distance of 976 at a posi- 
tion angle of 33° East of North. Additionally, component A has 
a binary c ompanion confirmed through lunar occultation mea- 
surements ( Dunha m! 1 977t iRi chichi et all 19991) . These measure- 
ments show evidence for a binary orbit that changes the compan- 
ion's relative position significantly over 20 years (offset 45mas, 
PA 300° in 1977, offset 14mas, PA 120° in 1999). 

The visual binary was resolved in TIMMI2 and VISIR N 
band observations at a separation of 9V 8 3 ±0705, PA 30°±1° 
(measured in VISIR image). The secondary component is 
also resolved at M (central wavelength 4.6/iin) with a flux 
of 119±35mJy. We do not resolve the separate components 
of visual component A. Calibration uncertainty was high 
in the VISIR observations, as indicated by the low flux 
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Fig. 2. Top left: The SED of HD 71155 with a sin- 
gle blackbody fit to the excess emission at a tempera- 
ture of 120K (symbols > 5^m are shown after subtrac- 
tion of the photospheric contribution). This suggests an 
offset of 34AU (0790) from the star. This fit is marked 
on the extension limits plot left. The shaded regions are 
where the disc(s) is (are) expected to lie. The limits sug- 
gest an alternative two-disc model which also agrees well 
with the IRAS 12 yum data, as shown top right, is more 
likely. Values of F<a sc are taken from the SED fits for each 
model, using scaling to the MIPS photometry and the un- 
certainties therein to give the errors. 



measured on the primary (Table 0, and thus we do not 
photometrically confirm the excess. The flux of the binary 
given in Table [2] and shown on the SED (Figure [3]) was 
scaled to the expected primary photometry (so multiplied by 

^primary, expected/ F primary, detected)- The visual binar y was 
fitted with a K7 spectral type to fit the JHK photometry from the 
2MASS catalogue. Adopting the parallax distance to the primary 
of 29±2 pc the luminosity of the binary is 0.12L o , consistent 
with a luminosity of 0.1L o typical for K8-type stars. 

After subtracting the primary and binary contributions pre- 
dicted from the SED fitting the IRAS measurements still indi- 
cate significant excess at 12 and 25/im (see Figure [3] and Table 
Q]). There are no nearby 2MASS or MSX sources likely to be re- 
sponsible for the excess measurements in the IRAS results, the 
source is not in the galactic plane (b = 13.2), and no additional 
sources are detected in our field of view, thus the IRAS excess 
is not likely to be due to a background source (limits on back- 
ground sources listed in Table 0. No evidence for extension is 
found in the images. We should have detected discs larger than 
6. 1 AU assuming a disc flux of ~260mJy at Q from a dust temper- 
ature fit of 420K (scaled to IRAS 12 and 24/um photometry, error 
on the blackbody fit is 80K from errors on excesses measured 
in IRAS). Blackbody grains at this temperature would be at an 
offset of 2.2AU, or if the grains are small and lie at 3x black- 
body offset (as discussed in sections 3.3 and 4.1) they would be 
at 6.6AU, just beyond the limits on disc extension. The excess 
emission and SED fit should be confirmed before we can inter- 
pret these limits in terms of constraints on the emitting grains. 

HD 23281: HD 23281 was first identifie d as a possible host 
of mid-infrared excess bv lShvlaia & Ashokl (120021) . IRAS pho- 



tometry at 25//m is indicative of excess emission at 3.8<x signif- 
icance (TableQ]). Our photometric results do not confirm the ex- 
cess on this target (Q band photometry is consistent with excess 
but at a < 3cr level, see Table0. No additional objects are seen in 
the field of view, with a limit on undetected sources of < 4 mJy at 
N. There are no bright 2MASS or MSX sources nearby that may 
have been caught in the IRAS beam and could be responsible 
for source confusion. There is no indication of extension found 
on any of the images of this object. Narrow discs at >9.5AU or 
broader discs >15.9AU should have been detected as extended 
emission in the VISIR Q band imaging (assuming f disc = 24mJy 
from SED fitting, see FigureO . We fit the IRAS 25^m photom- 
etry and VISIR upper limits on excess with a dust temperature of 
210^gjjK . Blackbody-like grains at 21 OK would lie at an offset 
of 5.4AU, consistent with the limits on extended emission. Very 
small grains that could lie at 3x this offset (see sections 3.3 and 
4.1) should have been detected as extended emission, if Fdisc= 
25mJy although the level of excess and dust temperature fit are 
currently too uncertain to allow constraints to be placed on the 
dust properties. 

HD 75416: HD 75416 (77 Cha) was identified a s a po ssi- 
ble mid-infrared excess host by Mannings & B arlowl (Il998l) in 
their study of the IRAS catalogues. It has significant excess at 
12/im (Table 0. Our observation (TIMMI2, N band) did not 
confirm the excess (detected 204+88mJy, predicted stellar fl ux 
is 144±6mJy). MIPS photometry at 24/um (R ieke et al.ll2005l) is 
in good agreement with the IRAS detection at 25 fim (MIPS 128 
± 13 mJy photos phere 34±2 mJy , IRAS 117 + 38 mJy photo- 
sphere 31±2mJv). ISu et al] ((2006) presented new 24 and 70/im 
photometry which also confirms the excess. The limit on any 
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Fig. 3. The observations of A Gem and HD 2328 1 . Top: The results of observations of A Gem. Bottom: The results of observations 
of HD 23281. Symbols >5/mi are shown after subtraction of the photosphere (not for the binary). Limits on extended emission 
obtained from the imaging are shown in the left-hand plots. The shaded region marks the area that a disc would not be detected as 
extended emission. The disc location marked on these plots is taken from the SED fit for disc flux, and from assuming that the dust 
is blackbody-like to determine a radius. 



background source within the TIMMI2 field-of-view is <48mJy 
(3cr at N), and thus it is very likely the excess is centred on the 
source. There is no evidence of extended emission in the image, 
and thus we can place a limit on the disc radius of <55AU on 
a disc of flux 55mJy in the N band (based on a fit to the MIPS 
and IRAS photometry). The fit to the excess at a temperature of 
250K would translate to 1 1 . 1 AU for blackbody grains, consistent 
with the non-detection of extension in the TIMMI2 image (error 
on blackbody temperature 60K from errors on excess emission). 
Even very small grains at 3x the blackbody radial offset (see 
section 4.1) would not be detected as extended emission. 

HD 141795: Th i s star was listed as an excess candidate 
bv lShvlaia & Ashokl d2002l) . The IRAS 25^m measurements of 
this source's photometry suggests an excess of 141 ± 28mJy. 
Calibration uncertainties prevent a photometric confirmation of 
the excess in the TIMMI2 and Michelle observations of the 
target (Table |2). Background and companion sources are ruled 
out at a level of 49mJy (N band, TIMMI2) and 8mJy (Q band, 
Michelle), and thus any excess is not likely due to detection of 
an additional source in the IRAS beam. 

The images show no evidence for extension, and we place 
limits on the extension of a disc with flux of 146mJy at Q (es- 
timated from the 25^m excess measurement and 12 and 60/im 
upper limits) of <6.2AU (Figure 0). This is consistent with the 
SED fit which uses a dust temperature of 250K putting the dust 
at 4.2AU for blackbody grains. However, the true dust temper- 
ature and disc flux at Q is highly uncertain (error on blackbody 



temperature fit is 70K). Photometric confirmation of the excess 
is necessary to confirm the limits provided by the extension test- 
ing. 

5. Resolvability of discs in the mid-infrared 

We now consider what mid-infrared debris discs could be re- 
solved with currently available instruments and future instru- 
mentation. We use the e xtension test i ng m ethod described in 
section 3.3 and in detail in lSmith et al.l d2008l) with PSF size and 
sensitivity appropriate to each instrument considered to deter- 
mine the limiting disc parameters for resolution (disc size and 
flux for different geometries). A short description of the param- 
eters used for each instrument considered are given in the sub- 
sections below, and the limits shown in Figure [5] 

On each plot we also show a representative sample of A star 
debris discs for comparison with the determi ned limits (overplot- 
ted with filled circles). This sample is from lWvatt et al.1 (|2P07b) 
and shows A star discs detected at 24 and 70/mi or 25 and 60yimi. 
Disc radii and flux levels are taken from fits to the excess emis- 
sion with a single temperature blackbody as described in that 
paper. There are several uncertainties inherent in this fitting. The 
radii determined for the discs assumes the emitting grains be- 
have like blackbodies, when in reality small grains which are in- 
efficient emitters may dominate the emission and the dust could 
be up to 3 times further from the star than this blackbody ra- 
dius. Multi-temperature fits to the excess are possible. Different 
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Wavelength, /xm Radius, arcsec 

Fig. 4. The results of analysis of the observations of HD 75416 (top) and HD 141795 (bottom). Left: The SEDs with symbols 
at >5jum shown after subtraction of the photospheric contribution at this wavelength. Right: The limits on disc location given by 
non-detection of extension in the images. The grey areas show possible disc locations for any of the geometries tried. The filled 
circle gives the predicted location from the SED fit. 
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38678 


<T Lep 


18 


0.43 


475 


0.18 


400 


Moerchen et al. (2007b) 


39060 


/JPic 


18 


1.59 


2688 


2.59 


4336 


Telesco et al. (2005) 


109573 


HR4796 


18 


0.51 


1307 


1.04 


807 


Telesco et al. (2000) 


141569 




18 


0.45 


549 


0.63 


623 


Fisher et al. (2000) 


181296 


zyTel 


18 


0.53 


142 


0.5 


62 


Smith et al. (2009a) 


39060 


/?Pic 


25 


1.24 


8005 


2.59 


6960 


Telesco et al. (2005) 


172167 


Vega 


24 


14.7 


1271 


11 


1500 


Su et al. (2005) 


216956 


Fomalhaut 


24 


8.7 


939 


20 


700 


Stapelfeldt et al. (2004) 



Table 4. The predicted and measured disc parameters of sources with resolved debris discs in the mid-infrared. 



temperatures could represent different grain populations at one 
radial location, or could indicate dust at seve ral radial locations 
(as is the case for rj Tel, Sm ith et al.ll2009al) . In such cases the 
predicted radial location and the level of disc flux for each com- 
ponent would be different from the simple single temperature 
fit shown here. The level of disc flux predicted at wavelengths 
other than 24 and 70/vm (or 25 and 60/im) may be incorrect even 
in the case that a single temperature is an accurate model for 
the emission. This is particularly true if spectral features are in- 
volved, for example several spectra of debris disc targets with 
IRS on Sp itzer have shown strong silicate features in the N ban d 
(see, e.g. lRieke et"aL1l2005b IChen et al.ll2006l: iLisse et al.ll2009l) . 
The effect of these uncertainties can be seen in the discs already 



resolved (shown in red in Figure|5]l. The resolved disc locations 
and fluxes are shown by asterisks and listed in Table|4] The value 
of /^obs/^pred (observed disc radius / predicted disc radius) is as 
high as 2.3 (for Fomalhaut) for the restricted set of discs resolved 
in the mid-infrared. This ratio is as low as 0.42 for HD38678, 
which may have a multiple component disc (as suggested by 
iFitzgerald et al.ll2007l) which was incorrectly fitted with a single 
temperature. As a final note of caution, the population shown on 
these plots is only a sample of known discs detected at 24 and 
70/mi and thus may not be truly representative of the population 
of discs at 10 and 18/im (e.g. hot discs detectable at 10/iin may 
not be detected strongly at 70/im). Discs around Sun-like stars, 
which will in general be smaller than the A star discs (as dust 
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must be closer to cooler stars to heat to mid-infrared tempera- 
tures) are also excluded from the sample shown. These plots can 
be used as a guide to the best sources to include in future ob- 
servational programmes aimed at resolving mid-infrared discs, 
but only through such resolution can the true disc parameters be 
known. 

5. 1 . Gemini instruments 

We consider the detection limits achieved in 2 hours of observ- 
ing at 18 and 25ji/m (2 hours on-source; after overheads and re- 
peated standard star observation to monitor the PSF total ob- 
serving time approximately 8.5 hours). The PSF model used 
is a Gaussian with FWHM of 076 at 18/mi (typical of 18/itn 
observations presen ted here) and (X'72 at 25/mi (taken from 
iTelesco et al.l 12005 b. The point source sensitivity follows from 
the detection levels found in the Of! 5 apertures at 18jt/m and is 
1.8mJy in 2 hours on source; extrapolation to 25/vm (sensitiv- 
ity 4.8mJy in 2 hours) assumes a factor of 8/3 brightness in- 
crease needed for a source to achieve the same signal-to-noise 
in the Qb filter (25//m) of TReCS as in the Qa filter (18/an), 
as outlined on the Gemini website. The detection limits for ex- 
tended disc emission were determined in the same way as the 
extension limits for the observations presented in this paper. 
Models of disc+star emission (with disc geometries and disc 
flux as described in section 3.3) convolved with the model PSF 
were treated as model images, and subjected to extension test- 
ing. Point-source subtracted model images were tested for sig- 
ni ficant residual emi ssion in optimal testing regions as described 
in Smith etal. (2008). Emission above 3cr significance was re- 
garded as a detection of extended disc structure. The resulting 
limits (Figure [5] top line) show that the best targets for resolved 
disc imaging campaigns are those that have already been re- 
solved (as for Figures 1-4, the region above the lines represents 
the disc parameter space that would result in a significant de- 
tection of extension acc ording to the met hod outlined in sec- 
tion 3.3 and in detail in Smi th et al.ll2008h . This plot was used 
to identify rj Tel, whose excess was independently confirmed 
in the TIMMI2 data presented here, as a prime tar get for 8m 
resolu tion. The resulting observations presented in Smit h et all 
(2009a) resolved the outer disc component and highlight the 
utility of this technique. Of the known A star debris discs pop- 
ulation few sources remain that could be resolved in reason- 
able observing times with current instruments at 18/mi (those 
most amenable to resolved imaging are HD 19356, HD 139006 
and HD102647 from current predictions of the disc parameters). 
More discs could be resolved at 25//m, although conditions suit- 
able for 25//m observing are more rare. 

5.2. MIRI on the JWST 

The James Webb Space Telescope is due to be launched in 
2013. MIRI, the mid-infrared instrument, will provide vastly 
greater sensitivity to debris discs in the 5-27/im spectral range. 
Since the primary mirror of the JWST consists of 18 hexago- 
nal segments with the combined collecting area of 25m 2 (sim- 
ilar to a 6m circular primary), the shape of the mirror seg- 
ments is reflected in the shape of the expected PSF. Thus the 
PSF was taken from the JWPSF tool for an A5 type star (tool 
at http://www.stsci.edu/jwst/software/, see Figure [6] for PSF). 
Noise from possible PSF uncertainty was calculated as the dif- 
ferences in residual flux emission found when using different 
PSF models as determined using the different OPD files avail- 



Fig. 6. A model 

i of the MIRI PSF 

at 18yum made 
using the JWPSF 
tool. The model 

^PSF is shown 
after subtraction 
of the best fit- 
ting Gaussian 
to show the low 
level structural 
features resulting 
-2 o 2 from the seg- 

mentation of the 
primary mirror. 



able (files represent modelled mirror misalignment and aberra- 
tions). We also tested PSF models based on AO type stars and 
found no significant difference in the limits achievable . The opti- 
mal re gions for testing for extension as described in Smit h et al.l 
(2008) and used to determine the extension limits on sources 
in this paper (see section 3.3) were modified to exclude pix- 
els close where differences in the PSF model from different 
OPD files were high (>50% the level of the signal of the PSF 
model in a pixel). The point source sensitivities were assumed to 
be 0.7juJy at lO/vm, 4.3/vJy at 18jum and 28yuJy at 25/mi (10cr 
in 10,000s) following the guidelines on the JWST webpages 
(http://www.stsci.edu/jwst/science/sensitivity/). These sensitivi- 
ties were converted to 1 hour's on-source integration by assum- 
ing that signal-to-noise varies as time 05 . Overheads have not 
been included in this calculation. 

The resulting limits on resolved disc parameter space are 
shown in Figure[5] middle panel and bottom left for face-on discs 
at lOjjm. The increase in sensitivity over current ground-based 
imaging is clear from the low levels of disc flux for which detec- 
tion of extended emission is possible in only 1 hours observing. 
Almost 100% of the A star discs detected at 24 and 70jum should 
be resolvable with MIRI, although the resolution of the discs 
close to the inner radius limit will strongly depend on accurate 
PSF calibration (and of course the caveats relating to disc flux 
/ radius predictions from SED fitting must be considered). The 
resolution of edge-on discs is a strong function of position angle 
as the PSF is not circularly symmetric. Discs which lie along the 
direction of the corners of the hexagonal shape (see Figure[6]l are 
more difficult to detect in residual emission as these regions see 
more noise resulting from mirror misalignment and aberrations; 
observing at different position angles could mitigate against this 
issue. The 'bump' in the detectability limits for face-on ring 
models also arises from the hexagonal feature of the PSF. Wider 
rings and discs lying edge-on are less strongly effected by this as 
their detection depends on less confined radial locations. A four 
quadrant phase-mask (4QPM) coronagraph will be offered at 3 
wavelengths; 10.65, 11.4 and 15.5/itn (w avelengths optimised 
for planet detection. lBoccaletti et al.ll2005h . We approximate the 
effect of including a 4QPM at lOyum (bottom panel, Figure[5} by 
increasing the sensitivity of the observations by a factor of 250 at 
Q3A/D and 50 at 5/l/D, falling to a factor of 1 at lOAfD (where D 
is 6.5m ). These values are based on Figure 8 of iBoccaletti et al.l 
(2005). Including these sensitivity improvements allows the de- 
tection of discs of flux down to a limit of ~9^Jy if the disc is at 
the optimal detection radius ( 1'.'44, minimum of dot-dashed line 
in Figure [5] bottom panel). Without the coronograph the mini- 
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Fig. 5. Predictions for the resolvability of discs with current 
and future instruments. See text for details of model limits and 
disc properties. Lines represent 3<x detection limits as described 
in section 3.3. Different source geometries are shown by the 
different lines and described in the legend (see also section 
3.3). Sources in red have already been resolved at this wave- 
length. Circles mark disc parameters estimated from SED fitting. 
Asterisks mark the true disc fluxes and radial offsets seen at this 
wavelength for the resolved sources. Errors from PSF variation 
for Gemini instruments and the E-ELT are approximated by us- 
ing PSF model with FWHM varying from to 8+d8 (see text for 
values of 6). This uncertainty is assumed to be at the 10% level 
(d6/8 = 0.1) consistent with our 8m imaging data. PSF variation 
is modelled by different simulated PSFs for MIRI. We also show 
the limits that could be acheived for a perfectly stable PSF in the 
top panel (compare solid and dotted lines). 



mum disc flux required for the detection of extended emission is 
0.01 3mJy. A Lyot mask optimised at 23/mi will also be provided, 
and will be used primarily for the detection of cold circumstellar 
discs. However, due to the large opaque mask of the Lyot object s 
at < l"cannot be detected with the Lyot dBoccaletti et alj |2005). 
As we can explore within this radius with our PSF subtraction 



method we do not include the effects of a coronagraph in our 
predictions at longer wavelengths (middle panel, Figure|5]l. 
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5.3. METIS on the E-ELT 

The European Extremely Large Telescope (E-ELT) is currently 
planned to have a 42m dish and to start operation in 2018. In the 
mid-infrared the current proposed first generation instrument is 
METIS, which will cover the 3-13^m ra nge (although form erly 
Q band was proposed for this instrument, Br andl et al.l2008l) . We 
scale our PSF models from 8m instrumentation to a 42m dish 
giving a FWHM of 0706 at lOfim. The sensitivity for imaging 
was assumed to be 8^/Jy in 1 hour on source (10cr) as taken from 
the METIS webpages (http://www.strw.leidenuniv.nl/metis/). To 
approximate the effect of the inclusion of a coronagraph we as- 
sumed that the sensitivity would be improved by a factor of 10 
at our hard limit of disc detectability (07006 where we adopted 
a PSF variation of O.lxFWHM which is typical of the level 
of variation we see in our 8m imaging data). This factor was 
taken from the increase in sensitivity predicted for planet de- 
tection when including a f our quadrant phase mask on top of 
PSF subtraction as given in lBrandl et al.1 (120081) . For larger radii 
we assumed the sensitivity improvement would decrease linearly 
with radius to a fixed level of 1 at >075 (~8 times the FWHM 
of the assumed PSF). This model is highly over-simplified and 
does not take into account many issues such as how sensitivity 
improvements may differ from point source predictions for ex- 
tended emission . We note that work currently being undertaken 
by Eric Pantin and collaborators to determine the sensitivity of 
METIS to disc and planet detection will adopt much more real- 
istic PSF and coronagraph models. 

The predictions for resolvable disc parameter space with 
METIS at lOjt/m are shown in Figure [5] bottom left. The com- 
parison with MIRI on the JWST is representative of the differ- 
ent strengths of ground-based and space-bourne instrumentation. 
Larger discs with lower surface brightness will be ideal targets 
for MIRI observations, but a large dish like the E-ELT will be 
needed to resolve very small discs (<073), or indeed structure 
within larger discs. Although the A star disc sample shown on 
this plot does not fill much of the small scale region, this is be- 
cause discs detected at 24 and 70/jm are cool and thus at large 
radii. The inner components of multiple dust population discs or 
dust at 10/um around Sun-like stars will live in these small spa- 
tial regions (see earlier discussion). This 10/jm disc population 
is poorly known as studies of these discs have been limited by 
calibration accuracy (specifically we cannot detect discs by clas- 
sical aperture photometry fainter than the level of accuracy with 
which we know the stellar emission at this wavelength, which 
is typically 10% of F*, although interferometric techniques can 
allow this limit to be surpassed). 

METIS on the E-ELT would enable the discovery of disc 
populations which cannot be detected photometrically (see 
above) through the detection of extended emission at 10/um. 
At the optimal detection radius (0705) METIS should be able 
to detect extended discs at the level of 6//Jy. We can compare 
this value to the flux expected for an exozodiacal cloud around 
nearby stars. If a star has 1 zodi of emission between 0-3AU 
with constant optical depth then it is most likely to be resolved 
in observations with METIS at 0705. We therefore calculate the 
detectability of exozodiacal emission by assuming this is equiva- 
lent to the detectability of a ring at the optimal rad ius. Taking the 
optica l depth of the zodiacal cloud, t = 5 x 10~ 8 (iDermott et alJ 
2002a), and assuming a ring of width dr = 0.5r gives a fractional 
luminosity of / = 0.5(dr/r)r = 1.25xl0 -8 . If we consider a sun- 
like star at lOpc, the optimal detection radius would be centered 
at r — 0.5AU. Adopting the blackbody temperature for this dust 
(T = 278.3/ -0% see section 4.1) then the observed flux from 




800 



Fig. 7. Plot of age versus fractional excess from SE P fitting at 
24 fim. The line shows a t~ l dependence, taken from Rieke et al. 
(2005), and this plot can be directly compared to their Figure 2. 



ex o-zodiacal du st should be 10 fi Jy at 10/-/m (using equation 6 
of Wyatt 2008: F v = 2.95 x IQ-™ B V (A, T)fr 2 /d 2 where B v is 
the Planck function, A is the wavelength of observation assumed 
here to be 10/um and d is the distance to the star here assumed 
to be lOpc) with higher flux for closer stars. Thus we could ex- 
pect to resolve discs down to <1 zodi out to lOpc. This result 
is consistent with the previous expectations of the performance 
of METIS (as discussed in section 2.3.3 of lBrandl et ail [2008, 
METIS is expected to be able to resolve the exozodiacal emis- 
sion in the 1AU region around stars at <10pc.) We will also be 
able to resolve details of the strucutre of the few bright discs 
already known at 10/rni, such as those around HD69830 and r\ 
Corvi (K0V and F2V) which a re believed to have dust in the 
terrestrial planet regions (~lAU. ISmith et al.ll2~009bl) . 

6. Discussion 

This sample contains 7 sources with excess infrared emission 
confirmed either in this paper or with Spitzer data (Rie ke et al.l 
120051: IChen et alJl2006h . and a further 3 sources with excess 
emission in IRAS requiring confirmation. The SED fitting in- 
dicates that these objects are surrounded by dust at a distance of 
between 2-60 AU (or alternatively two temperatures of dust at 2 
and 61 AU for HD71 155 and 4 and 24 AU for 77 Tel). These re- 
gions are those in which we might expect the formation of giant 
planets, and so it is important consider how the existence of this 
dust emitting in the mid-infrared adds to our current understand- 
ing ofdiistdi^ibtrfions in circumstellar regions. 

iRieke et aT] (120051) looked at a sample of 266 A-type stars 
between 5-626 Myr old with MIPS at 24 and 70/im and exam- 
ined the relationship between fractional excess and the age of the 
central star. They found that the upper limit of excess emission 
generally fell off as time~ l for the stars with detected excess. 
Assuming the fits to the SED profiles presented in this paper and 
plotting the predicted 24/um excess emission compared to stel- 
lar flux versus age it is clear that the results presented here are 
in-line with this relationship (Figure[7). Combining this with the 
age spread in sources suggests we have a representative sample 
of A star debris discs. 

HD71 155 and 77 Tel, the two systems with resolved dust pop- 
ulatio ns, have bo t h bee n found to have multiple disc compo- 
nents. |CheneLal. (2009) summarised the debris systems which 
have, through SED fitting or resolved imaging, been identified 
as multiple component discs. Adding HD71155 to this sample 
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brings the total number to 12 kn own "Solar System analogues" 
(as defined in Ch en et al.ll2009l systems with multiple compo- 
nent discs). The FEPS survey on Spitzer surveyed 328 solar- 
like stars finding ~10% have 70/mi emission indicating the pres- 
ence of cold debris, and of these 1/3 have SEDs that are bes t 
fit by multiple temperature excesses (Hillenbra nd et al.l 12008). 
From our sample of 10 debris disc sources (3 of which, A Gem, 
HD23281 and HD141795, still require photometric confirma- 
tion), 2 are multi-component discs (20%; consistent with the 
1/3 rate given our small sample size). These results are further 
evidence that multiple component or extended discs are com- 
mon. The traditional view of debris discs is that of a ring of 
planetesimals residing outside any planetary system producing 
cold dust, analogous to the Edgeworth-Kuiper belt (EKB) in the 
Solar System. Multiple component discs could be seen as analo- 
gous to the Solar System, although of course the extrasolar discs 
detected to date are much brighter than the asteroid belt and 
EKB (Lir/L* > 10~ 4 for detected extrasolar discs; asteroid and 
Edgeworth- Kuiper belts L™/L+ = 10~ 8 - 10' 7 and 10" 7 - 10~ 6 
respectively, iDermott et aTl2002bl IStern & Colwelill 19971) . 

For sources with large amounts of em ission from close to 
the star, the origin of the dust is uncertain. IWvatt et aD (l2007al) 
presented a model based on the collisional evolution models of 
iDominik & Decinl (120031) that predicted a maximum dust lumi- 
nosity dependent on age for a disc at a given radius. This max- 
imum brightness occurs when discs just reach collisional equi- 
librium, in which the size distribution of the planetesimals in a 
disc is fixed as mass is transferred down through a cascade of 
collisions to smaller and smaller sizes until they are removed by 
radiation pressure. Very massive discs process their mass very 
quickly and are therefore short-lived, whereas sparse discs take 
a long time to reach collisional equilibrium but are not very lu- 
minous. This model has been shown to accurately recreate the 
A star debris disc population observed with Spitzer (Riek e~etafl 
120051) under assumption s of a distribution of initial disc mass, 
radius, and current age (Wvatt et al. 2007b). The predicted max- 
imum luminosity for each disc (or disc component) in this study 
is giv en in Table [3] (f ma x = 1.2 x IQ^r 1 ^t^L from equation 
20 in [Wvatt et al.ll2007bl where r is the radius of the dust belt 
from our fits and f age is the age of the system as given in Table 
[1} . Uncertainties in the model parameters mean that only if 
/iR//max > 1000 do we take this as evidence of transient emis- 
sion (level of excess cannot be produced by a collisional cas- 
cade). Most of our targets have /iR// max < 1000, with A Gem and 
the inner disc of HD71 155 exceeding this value, and HD141795 
having /iR// max ~ 550. Of these possibly transient sources only 
HD71155 has confirmed excess emission (see Section 4.2), and 
confirmation of the other discs would be required before specu- 
lating on their origin. 

At an offset of 2AU around HD71155 the collisional life- 
time of bodies is short and thus we interpret the emission as ev- 
idence for a transient dust producing event. A recent massive 
collision in an otherwise quiescent disc could produce a short 
lived increase in excess, although as such events are likely to 
be rare (and the resulting dust grains have short lifetimes) the 
pr obability of witne ssing such an event is low (see discussion 
in lWvatt et"aT1l2007al) . At an age of 169Myr (Table [TJ, ongoing 
terre strial planet formation coul d be responsible for the emis- 
sion f Kenvon & Bro mley 2004b), with collisions between plan- 
etary embryos resulting in large amounts of dust production. 
The cooler belt at 61AU could represent a parent population 
of the hot dust emission as well as producing spatially coinci- 
dent dust dominating the excess at longer wavelength s (a similar 
possibility exists for the older Sun-like star r\ Corvi, Smit h et all 



2008. l2009bl) . However, the transport mechanism to get the dust 
to 2AU is unknown. A dynamical instability like that thought 
to have caused the Late Heavy Bombardment (LHB) in the 
Sol ar System (triggered by the migra tion of Jupiter/Saturn, see 
e.g. lGomes etal]|2005l;lLevison et alj |2008). during which large 
amounts of debris from the E KB was thrown into the inner Solar 
System, could be responsible. iBooth et alJ J2009) concluded that 
such an event occured around at most 12% of Sun-like stars. 
Around A stars, where the 10/im excess has been shown to ex- 
hibit a fall-off proportional to time like the 24 /j.m emission (al- 
beit with a longer decay time, ISu et al. the statistics are 
less clear. Resolving the outer disc component would allow fur- 
ther examination of a possible link between the two populations 
(for example, emission spread inwards towards the hot dust pop- 
ulation rather than confined to a narrow could be evidence for 
a link). A further possibility is that the emission arises from 
dust grains not produced in collisions but in the sublimation of 
a population of comets or one Super-comet. These possibilities 
were explored as the origin o f the h ot dust population around 
HD69830 by iBeichman etafl d2005l) . who concluded that the 
continuous generation of small grains by a population of comets 
would require too large a mass reservoir to be the likely origin 
of the dust. A single massive comet (Sedna-sized in the case of 
HD69830) could relea se small dust grains ove r a few Myr if cap- 
tured into a close orbit Beichman et a Li (120051) . This mechanism 
could be responsible for the 24//m dust population in more sys- 
tems. 

Dust emission that is from a transient event will necessarily 
show temporal variation. The difference between the IRAS pho- 
tometry on HD3003 and the measurements taken with MIPS and 
those presented in this paper (see section 4.1) could be evidence 
of such evolution. If real, this variance could be a reflection of 
the binary nature of the system, with the orientation of the sec- 
ondary as it proceeds on its orbit changing the overall illumina- 
tion of the system. A determination of the orbit of the binary 
will allow this possibility to be checked. Alternatively, if the 
temporal variance reflects changing levels of dust or changing 
dust distributions then a transient origin is more likely. Taking 
the assumed circumstellar radius of 4AU for the dust (see sec- 
tion 4. 1) /iR// max = 379, a high level but one at which we would 
not state the emission must be transient con clusively (see abov e 
and detailed discussion of uncertainties in IWvatt et al.ll200"7ah . 
However, if the dust is of transient origin, the dust location needs 
not be stable, and the constraint of >14.4AU for the binary semi- 
major axis need not hold. In this situation the determination of 
the orbit of the binary could again greatly improve our under- 
standing of this system. Alternatively resolving the dust distri- 
bution (as would be possible with METIS on the E-ELT, see 
Figure |5]l could also inform our models of the stability of the 
system, and thus the likely origin of the dust. 

Of the sources considered in this paper only r\ Tel and 
HD7 1 155 have had the location of the dust population confirmed 
by resolved imaging. For the remaining discs constraints have 
been placed on the dust location, but degeneracies in the SED fit- 
ting in particular (summarised in Section 5) mean that resolved 
imaging is required to determine the true dust distributions in 
these systems, and so constrain models for the dust origin par- 
ticularly where transience is inferred. The predictions in section 
5 can be used to target sources most likely to be resolved with 
currently available instruments, but as shown in Figure [5] most 
currently known disc targets will require the use of MIRI to 
detect faint levels of extended emission or the high resolution 
of METIS on the E-ELT to resolve emission on small scales. 
High resolution will also be important for the detection of sub- 
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structure in the discs which could indicate the presence of plan- 
ets which will be important for ascertaining the nature of these 
systems and distinguishing between models for the origin of the 
dust. Evidence for planets in the dust distribution of debris disc 
systems include: sharp disc edges (as seen arou nd Fomalhaut, 
iKalas et al.l 120051) : clumps (Vega. IWvattl 12003: similar struc- 
ture may be observable in the EKB because of Neptune's res- 
onant Plutino popula tion if the disc were brighter); warps QS Pic, 
Augereau etal]|2001h : and asymmetries (HR4796A, Wya tt et all 
1999). 

7. Summary 

In this paper we have presented new observations of 1 1 early- 
type stars which have been proposed to be debris disc hosts 
based on their IRAS photometry. We have used TIMMI2, VISIR, 
Michelle and TReCS data to confirm excess emission and/or 
place constraints on debris discs for the observed sample. Our 
results are: 

- For HD 3003, HD80950, and 77 Tel our photometry yields 
an independent confirmation of excess emission around the 
target. Subsequent analysis of the HD3003 system indicates 
that if the dust lies in a stable region it must be circumstellar 
and the binary must orbit at a semi-major axis of >14.4AU 
(assuming blackbody grains). 

- Our data on 5 targets allow us to place quantitative limits on 
the location and level of emission of any dust in the system. 
For HD7 1 155 these limits allow us to determine that the disc 
must have multiple components. 

- We use simple disc models to determine the region of disc 
flux versus radius parameter space for which discs can be re- 
solved with currently available 8m mid-IR instruments. This 
technique successfully identified rj Tel as a resolvable disc, 
which was confirmed with TReCS (Smith et al. 2009a). 

- We predict the parameter space of resolvable discs that could 
be opened by future instruments MIRI on the JWST and 
METIS on the E-ELT. Spatially extended disc structures will 
be best observed with MIRI because of their lower levels of 
surface brightness, whereas discs close to their central star 
(within ~073), or those with structure on small spatial scales, 
will be prime targets for E-ELT imaging which would be able 
to detect emission below 1 zodi out to lOpc. 
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Appendix A: Additional sources 

A. 1. Observed sources with no new limits 

HP 31295: HP 312 95 was first identified by 
Sada kane & Nishidal (0~986) in their sample of Vega-excess 
stars as a star with infrared excess identified in the IRAS PSC 
(significant at 25/mi, see Table[TJ. TIMMI2 N band observations 
do no t confirm the excess (Table [2). R esults of Spitzer observa- 
tions (lJuraet al.ll2004t ISu et al.ll2Q06h place a limit of <70mJy 
excess at 8.5/xm, and the IRS spectrum shows no evidence of 
excess at <20 /mi. Further results presented in ISu et al, 1 (120061) 
showed that in MIPS 24 and 70/mi photometry the excess was 
confirmed. We find no evidence of background/companion 
sources in the TIMMI2 field-of-view (such objects limited to 
<33mJy). 

We fit the confirmed excesses with a blackbody at 80K (sim- 
ilar to the Chen et al. fit of 90K, our errors are + 12K from errors 
on the excess detections), which translates to a radial offset of 
52.6AU (1'.'4, see Figure |A~T1 for SED). The image shows no 
evidence for extension. The disc flux expected in the observed 
TIMMI2 filter from the SED fit is <lmJy; longer wavelength 
observations would be required to resolve/limi t the d isc exten- 
sion around this source. IMartmez-Galarza et alj (l2009h explored 
the possibility that this A Bootis star could be interacting with 
the ISM, and found that such interaction could produce the ob- 
served excess. With the current data it is not possible to dis- 
tinguish between this theory and that of a circumstellar disc, 
although the fact that this star lies within the local bubble (at 
37pc, see Table [1} means that the probability of the star lying 
withi n a cloud which could prod uce the observed emission is 
low (IMartmez-Galarza et al1l2009h . 

HD 38206: HP 38206 was first iden t ified a s a host of mid- 
infrared excess by Mannings & Barlow (1998) in their analy- 
sis of the IRAS catalogues. The TIMMI2 photometry presented 
here does not confirm t he excess (Tab le Recent MIPS ob- 
servations of this star dRieke et al.ll2005l) have confirmed the 24 
/mi excess, with flux of 115 + 12 mJy (expected photospheric 
flux 33+3 m Jy, good agreem ent with the IRAS measurements, 
see Table Q}- Su etalJ (l2006h list the MIPS 24/zm photometry as 




0.1 1.0 10.0 100.0 1000.0 

Wavelength, /xm 

Fig. A.2. The SEP of HP 23432 with symbols at > 5/mi show- 
ing emission after subtraction of photospheric contribution. The 
excess is most likely due to the reflection nebulosity in this re- 
gion. The TIMMI2 data points are highly untrustworthy due to 
the poor photometric conditions. Blackbody fits for illustration 
purposes are at 200K, 80K and 28K. 

107mJy at 24/mi and 342mJy at 70/mi (errors of 1.58mJy and 
12.87mJy respectively do not include calibration errors which 
are less than 5% at 24//m and 10% at 70/mi). The excess emis- 
sion is fitted with a blackbody at 90±10K (translating to a ra- 
dial offset of 48.4 AU, 0770, see Figure \KJ\ right and Table©. 
This fit suggests the disc flux at the wavelength observed with 
TIMMI2 is <lmJy; as for HP 31295, only longer wavelength 
high resolution imaging will be able to constrain or potentially 
resolve this disc's location and geometry. 

A.2. Not a debris disc candidate 

HP 23432: HP 2 3432 (asterope) was identified by 
lOudmaijer et al. as being amongst a sample of SAO stars 

with IRAS infrared excess. This star has an excess of 256 ± 43 
mJy at 12 /mi and 1159 + 42 mJy at 25 /mi (after subtraction 
of the photosphere). In addition it also has excess at longer 
wavelengths: 6533 ± 780 mJy at 60 fim; and 13078 ± 4979 mJy 
at 100/mi. 

The excess emission is not confirmed in the TIMMI2 obser- 
vations of HP 23432, as a flux of 141 + 62 mJy at 11.6 //m is 
found compared to an expected stellar flux of 136+7 mJy from a 
Kurucz profile fit (see Table©. The TIMMI2 data points plotted 
on the SEP (Figure IA.2b are shown with the calibration lim- 
its taken from the standards immediately before and after the 
science observation. The overall photometric errors are much 
higher, with a change of calibration factor over the course of the 
night of around 30%. Optical observations of this Pleiades mem- 
ber show i t to lie close to a diffuse reflection nebula Ced 19h 
ICederbladl (1 19461) . The shape of the excess emission spectrum 
(fit here with blackbody emission at 200+30K, 80+25K, and 
28±3K, where errors on the temperatures arise from fitting the 
excess emission within 3<x) suggests that this is the result of the 
emission from the reflection nebul a and no t from a dust popula- 
tion centred on the star itself. iGorlova et ail(l2006l) found that the 
Spitzer 24/im observations of this source were contaminated by 
flux from the reflection nebula. Similar interactions with inter- 
stellar dust have been sh own to be responsible for excess emi s- 
sion towards other stars (Kal as et aDl2002l : iGaspar et aD l2008). 
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Fig. A.l. The SED fits of the excess emission of two stars with excess at longer wavelengths than that of the observations presented 
here. Left: HD 3 1295 with excess fit by a blackbody at 80K. Right HD 38206 with excess fit by a blackbody at 90K. 



